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ABSTRACT 

We investigate the luminosity evolution of the host galaxies of radio-loud AGN 
through Hubble Space Telescope imaging of 72 BL Lac objects, including new 
STIS imaging of nine z > 0.6 BL Lacs. With their intrinsically low accretion rates 
and their strongly beamed jets, BL Lacs provide a unique opportunity to probe 
host galaxy evolution independent of the biases and ambiguities implicit in quasar 
studies. We find that the host galaxies of BL Lacs evolve strongly, consistent with 
passive evolution from a period of active star formation in the range 0.5 < z < 
2.5, and inconsistent with either passive evolution from a high formation redshift 
or a non-evolving population. This evolution is broadly consistent with that 
observed in the hosts of other radio-loud AGN, and inconsistent with the flatter 
luminosity evolution of quiescent early types and radio-quiet hosts. This indicates 
that active star formation, and hence galaxy interactions, are associated with the 
formation for radio-loud AGN, and that these host galaxies preferentially accrete 
less material after their formation epoch than galaxies without powerful radio 
jets. We discuss possible explanations for the link between merger history and 
the incidence of a radio jet. 

Subject headings: galaxies: active — BL Lacertae objects: general — quasars: 
general — galaxies: evolution — galaxies: jets — black hole physics 



1. Introduction 



From Hubble Space Telescope (HST) studies of local galaxies, it appears that nearly 
all contain supermassive black holes (SMBHs) (Joseph et al. 2001; Bower et al. 2001; Sarzi 



- 2 - 



et al. 2001; Barth et al. 2001; Gebhardt et al. 2003). Combined with the predictions from 
hierarchical galaxy formation, this suggests that the formation and evolution of galaxies and 
SMBHs might be closely linked (Silk & Rees 1998; Haehnelt & Rees 1993; Kauffmann & 
Haehnelt 2000; Franceschini et al. 1999; Merrifield et al. 2000; Wang, Biermann & Wandel 
2000; di Matteo et al 2003). Circumstantial evidence also points to a direct link between 
the evolution of normal galaxies with the AGN phenomenon — the co-moving densities of 
both rise rapidly from today back to redshift ~ 2 (Schmidt & Green 1983; Boyle, Shanks 
& Peterson 1988; Shaver et al. 1996; Lilly et al. 1995; Driver et al. 1995; Cowie et al. 1996; 
Madau et al. 1996; Connolly et al. 1997; Cowie et al 2003; Hasinger 2003). Certainly it is 
plausible that AGN luminosity affects host galaxy properties, particularly star formation 
rates (e.g. de Young 1989; Rees 1989; Whittle 1992; Gonzalez Delgado 1995; Barthel 2001; 
Zirm, Dickinson & Dey 2003), and a number of studies have found a dependence of maximum 
AGN power on host galaxy luminosity (Smith et al. 1986; Vefon-Cetty & Woltjer 1990; 
Hutchings et al. 1984; Hooper, Impey & Foltz 1997; McLeod, Rieke & Storrie-Lombardi 
1999; Schade, Boyle & Letawsky 2000; Dunlop et al. 2003). 

It may be that most galaxies pass through one or more active phases at some point 
in their development (Cavaliere & Padovani 1989). If so, then comparing the luminosity 
evolution of AGN host galaxies to that of normal galaxies allows us to study the link between 
AGN activity and galaxy evolution. If galaxy interactions are an important trigger of AGN 
activity (e.g. Hutchings & Campbell 1983; Heckman et al. 1984; Yee 1987; Heisler 1991; 
Hutchings & Neff 1992; Canalizo & Stockton 2001) then there is a link between AGN activity 
and the assembly of elliptical galaxies, which should be reflected in the luminosity evolution 
of AGN host galaxies (Franceschini et al. 1999). 

Work on host galaxy evolution to date has largely been restricted to high-power sources 
due to the difficulty in finding low-power AGN at sufficiently high redshift for the evolution 
of galaxies to become significant. The host galaxies of powerful radio- loud AGN at high 
redshift appear to be brighter than general radio-loud host galaxies observed at low redshift, 
consistent with passively evolving stellar populations (Aragon-Salamanca et al. 1993; Lilly 
& Longair 1982; Kukula et al. 2001; McLure & Dunlop 2000; Lehnert et al. 1999). The same 
is not observed of brightest cluster galaxies (BCGs), which exhibit flat luminosity evolution 
over the redshift range to 1 (Aragon-Salamanca et al. 1993), and so are likely to be gaining 
in luminous mass. This is also tentatively the case for other early-type galaxies (Stanford 
et al. 1998; Lilly et al. 1995; Kauffmann et al. 1996; Bell et al. 2004). The host galaxies of 
radio-quiet AGN seem to exhibit similar flat or negative luminosity evolution, and appear 
to be gaining mass between redshift z ~ 2 and the present (Kukula et al. 2001; Rix et al. 
2001; Ridgway et al. 2001). 
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However, various biases make it impossible to confidently ascribe the observed evolution 
to the host galaxies themselves. First, it appears that maximum AGN power increases with 
host mass (Smith et al. 1986; Vefon-Cetty & Woltjer 1990; Hutchings et al. 1984; Hooper, 
Impey & Foltz 1997; Schade, Boyle & Letawsky 2000; McLeod, Rieke & Storrie-Lombardi 
1999; Scarpa & Urry 2001; O'Dowd, Urry k Scarpa 2002; Dunlop et al. 2003). This 
may induce apparent redshift evolution in flux-limited samples. Second, in powerful AGN, 
scattered light from the nucleus and extended, ionized gas contaminates host galaxy light. 
In radio galaxies, nuclear emission may be mistaken for host galaxy light (Dunlop & Peacock 
1993). This contamination increases with AGN power, which may again produce a spurious 
host luminosity-redshift trend. 

BL Lac objects are intrinsically low-power AGN, with much lower accretion rates than 
quasars (O'Dowd, Urry & Scarpa 2002) or powerful radio galaxies, yet the relativistic beam- 
ing of their jet emission means that they can be detected in significant numbers to z ~ 1, at 
which point galaxy luminosity evolution becomes measurable. Fortuitously, BL Lac objects 
overcome many of the potential biases outlined above. 

First, the apparent luminosity of BL Lac objects is dominated by magnified emission 
from their relativistically beamed aligned jets (Urry & Padovani 1995). This magnification is 
in turn dominated by the angle of the jet to the line of sight. Variation in jet angle washes out 
any trend between the power of the nucleus and the luminosity of the host galaxy (supported 
by absence of any correlation between host luminosity and beamed nuclear luminosity; Urry 
et al. 2000.) Second, the host galaxies of BL Lac objects are less likely to suffer from 
contamination by scattered light from their low-power nuclei, nor do they tend to exhibit 
the luminous, extended, ionized gas observed in more powerful AGN. 

Due to the (1 + z) A surface brightness dimming of the hosts (compared to (1 + z) 3 
dimming of the nucleus), the Hubble Space Telescope (HST) is needed to reliably determine 
BL Lac host galaxy parameters beyond z ~ 0.6. In this paper, we present new HST 1 
STIS images of high-redshift BL Lac objects. Combined with earlier z < 0.6 HST WFPC2 
observations, we measure the evolution of radio-loud AGN host galaxies, free from many of 
the biases implicit in previous studies. The observations and data analysis are described 
in §2. Results on the host galaxy properties are in §3, analysis and results for the host 
galaxy luminosity evolution are presented in §4, and the discussion of these results in §5. 
Conclusions are summarised in §6. 



1 Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc. under 
NASA contract No. NAS5-26555. 
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2. Observations and Data Analysis 
2.1. The Sample 

The new STIS observations are high-redshift BL Lac objects selected from the cycle 6 
HST snapshot survey of Urry et al. (2000). The original 132 proposed snapshot targets were 
chosen from six flux-limited samples. These were essentially all known BL Lac objects in 
complete samples in 1995. The 110 BL Lacs finally imaged in the snapshot program were 
chosen at random from the larger sample, and are an unbiased subset. The sample chosen 
for STIS imaging comprised the 9 least optically luminous BL Lacs with z > 0.6, spanning 
a range at least to z ~ 1. As observed nuclear brightness is dominated by variation in jet 
angle, this luminosity requirement does not introduce significant bias, but does dramatically 
increase our chances of detecting the host galaxy. Table 1 gives the final target list and 
exposure information, with the magnitudes and host galaxy limits derived from the WFPC2 
snapshot images. 

These targets were imaged with the STIS camera aboard HST, using the F28x50LP 
longpass filter. This configuration maximizes sensitivity to the host galaxy, sampling from 
5,500A to 10,000A. The absence in BL Lac objects of a big blue bump or extended line 
emission, as seen in quasars, makes it possible to use such a broad filter. 

2.2. Data Reduction and Calibration 

Four to five exposures were taken for each target, dithered in a square with edge length 
5 pixels. In each case, one exposure was short enough to ensure that the nucleus was not 
saturated. The data were reduced using the STIS package in IRAF's STSDAS suite. The 
bias and dark frames were subtracted and the images flat-fielded using the routine STIS 
reduction process. CRREJ was used to perform cosmic ray rejection. 

Most of the long exposures resulted in over-exposure and bleeding of the BL Lac nuclei. 
In each of these cases the affected pixels in the combined image were flagged and replaced 
with the pixels in the shorter exposure, normalised according to their exposure times. 

2.3. Sky Subtraction 

Accurate determination of the background level and its uncertainty is critical in this 
analysis, as an incorrect measurement can easily lead to the false detection of host galaxies. 
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The background was taken to be the median in the region of each combined imaged which 
was unaffected by any visible source. Such sources were masked well beyond their visible 
extent for determination of the background. 

The error in the background value was calculated from two contributions: the pixel- 
to- pixel Poissonian shot noise (<Jp isson), and the large-scale variation in the background 
level {(J i ar g e - sca i e ). The latter was taken to be the maximum of: standard deviation in the 
medians of a grid of 20x20 pixel squares, and the standard deviation in the medians of the 
four quadrants of the image. The final error in the subtracted sky value was then: 



Figure 1 shows the central 300x300 pixels of the combined, sky-subtracted images. Also 
shown are the contour plots of the regions surrounding the targets, with north indicated by 
the arrow head and east indicated by the arrow tail. 



Perhaps the most critical step in extracting the properties of the host galaxies is the 
accurate modelling of the Point Spread Function (PSF). Unlike ground-based PSFs, the HST 
PSF is relatively stable with time, and so synchronous PSF observations are not essential. 
However, the PSF does vary significantly across the surface of the chip, and this must be 
taken into account. In all images the target was within 20 pixels of the centre, so such spatial 
variations are minimized. 

Due to deficiencies in the Tiny Tim model of the STIS PSF (it fails to account for 
extended, scattered light, the STIS 'ghost loop' and an asymmetric bulge near the PSF's 
centre), we elected to model the PSF with archival images of well-centred stars taken with 
STIS in the F28x50LP filter. Highly over-exposed stellar images were used to model the 
PSF wings, while unsaturated images were used to model the central regions. Figure 2 shows 
the final composite stellar PSF. The 'ghost loop', a reflection feature in the STIS PSF, was 
found to vary significantly even with small positional offsets, and so we simply mask this 
feature in all of the analysis to follow. 

Given the potential inaccuracies in the STIS PSF extreme care was taken to account 
for its uncertainties. Deviations in the PSF arise from spatial variations, time variations, 
and spectral variations. The stars that were averaged to make our composite PSF have 
different spectral types, were taken at different times, and are offset from each other spatially 




2.4. Modelling the STIS Point Spread Function 
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(although all are still relatively close to the centre). Thus, the deviations observed among 
these PSFs should reflect the potential difference between the composite PSF and our data. 
We define the uncertainty in each pixel of the composite PSF model to be the maximum 
difference between that pixel value and the corresponding pixels in the component stellar 
images. This gave a conservative uncertainty map which was used in the two-dimensional 
host galaxy fitting. 

2.5. Host Galaxy Fitting 

Host galaxy properties were measured using two-dimensional model fitting. The models 
comprised the stellar PSF (§.2.4) centred on an analytic galaxy model. The host galaxies 
were modelled with the following profile types: a de Vaucouleurs r 1 / 4 law, an exponential disk 
r _1 law, and a composite de Vaucouleurs plus exponential disk model. The free parameters 
fitted for each host galaxy profile type were: PSF normalization, host galaxy normalization 
(two normalizations for the composite model), host galaxy characteristic radius (again, two 
parameters for the composite model), ellipticity, and position angle. 

Certain parameters were limited or fixed, as the quality of the data did not justify more 
detailed fitting. The point source was fixed at the centre of the galaxy, as this is observed to 
be the case all in low redshift BL Lacs (Urry et al. 2000). The galaxy models were limited 
to the three profile types described, as we are primarily interested in distinguishing between 
disk- and bulge-dominated morphologies. Modelling more detailed profiles (such as via the 
Sersic index) is beyond the scope of these data. 

Before fitting, images were carefully masked to eliminate all irregular features such as 
close companions and tidal structure, along with several pixels beyond the visible edges of 
these features. The centre-most few pixels of the nucleus are highly sensitive to sampling 
errors. The values of these pixels will depend strongly on exactly where in the pixel the 
point-like nucleus falls. If the model PSF is offset even slightly from this (by as much as 
0.05 of a pixel), large errors will occur in the fit, as these central pixels strongly weight our 
fitting statistic. To remedy this, we sum the pixel values within a central circle with a 7-pixel 
diameter, and treat this as a single data point. As a side-benefit, this process greatly reduces 
the potential error due to an inaccurately modelled PSF. 

In all cases but that of 2240-260, the fits indicated the presence host galaxies. The 
residuals of 0820+225, 1144-379, 1249+174, 1422+580 and 1533+535 were smooth and 
either circularly or elliptically symmetric about the nucleus. The residuals of 0138-097, 
0235+164 and 1308+326 were less convincing, so great care was taken to determine the 
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confidence of the detections. 



2.5.1. Detection Criteria and Upper Limits 

To determine whether the apparent host galaxy detections were significant, we fitted 
each image with the PSF alone, and compared the x 2 values of the best PSF-only fit to the 
best PSF + host galaxy fit via the F-test. We required that the model including the host 
galaxy be preferred over the PSF-only fit with 99% confidence. 

To make certain that we did not mistake improperly subtracted background light for 
the host galaxy, we also perform the fits with the background over-subtracted by 2a S k y (see 
§2.3), and required that the model including the host galaxy was preferred over the PSF-only 
fit with 95% confidence in this case. 

To ensure that none of the detected host galaxies were simply artifacts of a poorly- 
modelled PSF, we redid the fits using each of the individual stellar images as our PSF model 
instead of the combined stellar PSF, and required again that the host galaxy detected with 
each fit was preferred over the PSF-only fit with 99% confidence. For each alternate PSF, 
we again require the host galaxy be preferred at 95% with the background over-subtracted. 

After applying these requirements, it was found that five of the nine host galaxy were de- 
tected with high confidence. These were the host galaxies of 0820+225, 1144-379, 1249+174, 
1422+580 and 1533+535. Of the remainder, 0138-097, 0235+164 and 1308+326 yielded 
host galaxies in the initial fit, but these did not meet our confidence requirements when the 
background was over-subtracted and/or the alternate PSFs used. 2240-260 did not show a 
detectable host galaxy in any of the fits. 

Figure 3 shows the one-dimensional, azimuthally averaged profiles of both the data and 
the best fits. Galaxy profiles are presented only in cases where the PSF + host galaxy 
model is preferred over the PSF-only model according to the initial resolution criterion; ie. 
for 0138-097 and 0235+164 and the exponential disk fit of 1308+326, where the model was 
rejected only after increasing the sky background or after fitting with an alternate PSF, we 
still present the best-fit PSF + galaxy profiles. 

Figure 4 shows the images with PSF subtracted. For the detections we use the best 
PSF normalization from the de Vaucouleurs fit, and for the non-detections we use the best 
PSF-only fit. For the detections we also show the images with best-fit PSF + de Vaucouleurs 
model subtracted. 

As a last check, we take the ratio of light in a aperture of radius 3.5 pixels around the 
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centre to the light an annulus from 15 to 50 pixels, and compare to the same ratio for the 
PSF. In all five of the confident detections, the BL Lac image showed a significant excess 
of light in the annulus, indicating the presence of host galaxies. In the non-detections the 
excess was marginal. 

For the non-detections, we determined the upper limit on the host galaxy magnitude 
by finding the maximum amount of host galaxy which could be added before the PSF + 
host galaxy fit became worse than the PSF-only fit at the 99% confidence level by an F-test. 
In these cases, the galaxy model used was a de Vaucouleurs model with zero ellipticity and 
an effective radius of 10 kpc — slightly higher than the median found for the low redshift 
BL Lac hosts of 8.5 kpc (Urry et al. 2000). This is conservative, as it will give a slightly 
higher upper limit than for the median r e . 

2.5.2. Uncertainty Analysis and Resolving Morphologies 

The uncertainty in each fitted parameter was determined by mapping x 2 around the best 
fit. To account for the uncertainty in the background subtraction, the fits were performed, 
and the x 2 volume was mapped with the background both under- and over-subtracted by 
cr sky (see §2.3). The uncertainty in each parameter was then the maximum displacement 
from that parameter on the one-sigma hypersurface of the x 2 volume. 

For a statistically significant preference of one profile type over the other, we again 
employ the F-test, requiring that their minimum x 2 values differ with 99% confidence. In 
most cases, the de Vaucouleurs model was marginally preferred over the exponential disk 
model, however in only two of these cases was the preference significant — for 0820+225 and 
1422+580. In no cases was the exponential disk profile preferred over the de Vaucouleurs 
profile with high confidence. 

As might be expected, the composite de Vaucouleurs + exponential disk model gave a 
better fit than either the de Vaucouleurs or exponential disk model alone. Typically, the 
improvement in the fit was marginal — at less than 90% confidence according to the F-test in 
all cases but 0820+225, for which the three-component model was preferred over the best-fit 
two-component model with ~90% confidence. 

A number of the de Vaucouleurs model fits gave unphysically large effective radii, sug- 
gesting that the fit may be trying to compensate for an extended, disky component. Restrict- 
ing the PSF + de Vaucouleurs model to more physical scale lengths — requiring r e < 20 kpc 
- then the composite galaxy model is preferred over the new best-fit PSF + de Vaucouleurs 
model for 0820+225 at the 99% confidence level. If we apply the same condition to the rest of 
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the sample, both 1422+580 and 1533+535 prefer the composite model at ~ 95% confidence, 
and 1144-379 at the ~ 90% level. 

Figure 3 shows x 2 contours projected onto the mhost^e plane for the de Vaucouleurs 
and exponential fits of the resolved sources. 

2.6. K Corrections and Filter Conversions 

Cousin's R band absolute magnitudes are presented for purposes of comparison with the 
low redshift BL Lacs in Urry et al. (2000). These corrections were performed by assuming 
SEDs for both the host galaxies and the nuclei. The models of Bruzual & Chariot (1993) 
were used for the host galaxy SEDs, calibrated to the galaxy colours reported in Fukugita, 
Shimasaku & Ichikawa (1995): a Hubble-type Ea spectrum (B — R = 1.57) for correcting 
the de Vaucouleurs fits, and a Hubble-type Sab (B — R = 1.34) for the exponential fits. A 
power law SED with a = — 1 was assumed for the BL Lac nuclei. 

For the filter conversions, the redshifted SED was convolved with the transmission 
curve for the F28x50LP filter and then normalised to the best-fit host galaxy magnitudes. 
K corrections were performed assuming H = 70 km s _1 Mpc _1 , VLm = 0.3, Vt\ = 0.7 (the 
cosmology used throughout). 

The level of uncertainty in the host galaxy spectral corrections is potentially very 
large. For sources at these redshifts, the F28x50LP filter samples well below the rest- 
frame 4000A break, and hence is highly sensitive to the star formation history of the galaxy. 
Uncertainties may be as large as 0.6 mag if the high-redshift host galaxies have colours closer 
to those of spiral galaxies, rather than the red colours observed at low redshift (Scarpa et al. 
2000a), or even as large as 1 mag if they contain a strong starburst component. The quoted 
absolute magnitudes should be interpreted with this in mind. 

Galactic extinction was corrected using the relation to HI column density from Shull & 
Van Steenberg (1995): E(B - V) = logN (HI)/ (21.72 cm^mag- 1 ) 

3. Host Galaxy Results 

The results of the fits for the single-component galaxy models and the PSF-only model 
are given in Table 2. Table 3 shows the results for the composite galaxy models. 
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3.1. Host Galaxy Brightnesses 

The average absolute R-band magnitude of the five resolved host galaxies — 0820+225, 
1144-379, 1249+174, 1422+580, and 1533+535 — for the de Vaucouleurs model fit, K cor- 
rected assuming an early-type spectrum, is Mr = —24.3 mag. This is much brighter than 
the average magnitude of M R = —22.85 mag measured for the lower redshift BL Lac host 
galaxies in Urry et al. (2000). This strongly suggests that there has been significant lu- 
minosity evolution between the two samples, even taking into account the presence of the 
(high) upper limits in this sample. The average magnitude remains high even if a much 
bluer spectrum is used for K corrections: —23.7 mag for a Hubble-type Sab spectrum, and 
—23.51 mag for a 1 Gyr-old burst population. 

3.2. Host Galaxy Morphologies 

We were able to resolve complex morphologies in two cases. Both 0820+225 and 
1422+580 clearly have bulge-dominated host galaxies. The effective radii derived from al- 
most all of the de Vaucouleurs fits are very large, although the sizable errors on these radii 
encompass more sensible, smaller values. Nonetheless, the derived radii of 0820+225 and 
1533+535 are unphysically large. These large scale sizes may indicate that some of the 
host galaxies in this sample may have substantial disk components as well as bulges. The 
bulge + disk host galaxy model is preferred with high significance for 0820+225, and with 
marginal significance for 1422+580 and 1533+535. 

Two of the unresolved sources, 0235+164 and 1308+326, show a significant amount of 
extended light which is fit well by an exponential disk component. This light can also be 
accounted for as resulting from an imperfectly subtracted sky, at the limit of the errors in 
the sky subtraction (and, in the case of 0235+164, at the limit of the potential variation 
in the PSF). As such, the two are considered non-detections. However, many of the host 
galaxies show the same excess of extended light above the best-fit de Vaucouleurs profile. It 
is possible that these 'non-detections' are actually disky hosts on large scales. These sources 
would require longer integration times to confidently resolve the host galaxies. 

3.3. Close Environments and Evidence of Interactions 

Most of the BL Lacs in this sample appear to have close companions, nearby irregular 
structure, or both. 1533+535 is the only source that has no such observable features. It 
is impossible to judge in the case of 0235+164 because of the known intervening system. 
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Discounting this BL Lac, five out of the remaining eight have apparent companions within 
a projected distance of 50 kpc, and four of the eight show signs of recent interactions, either 
companions within 10 kpc or irregular structure. 

This sample differs from the broader BL Lac sample of Urry et al. (2000), which ex- 
hibited close companions within 50 kpc in 47% of cases, rather than the 62.5% as observed 
here. Very few of the broader sample showed signs of recent interactions, compared to 50% 
of this sample. 

3.4. Notes on Individual Objects 

0138-097 The host galaxy of this BL Lac appeared to be marginally resolved using the 
composite stellar PSF. With two of the alternate stellar PSFs, however, it produced fits which 
were not sufficiently better than the PSF-only fit. Thus, it is considered unresolved. The 
difficulty in resolving this host galaxy is compounded by the presence of a bright companion 
object 1'.'44 (10.5 kpc) to the southwest of the nucleus, and more structure within 5". 

0235+164-' The field of this object is complex, with bright structure extending to ~ 3" 
(~ 30 kpc) from the nucleus. This structure has been well-studied in the past (e.g. Stickel et 
al. 1988; Yanni, York & Gallagher 1989; Abraham et al. 1993; Falomo 1996), and is known 
to be an intervening system with redshift z = 0.524. These STIS data resolve the complex 
structure in significantly more detail than previous observations, but we leave its analysis to 
a later work. 

Fortunately the extended structure is concentrated to the south of the nucleus, leaving 
most of northern half clear for our host galaxy analysis. A host galaxy appeared resolved 
using the composite stellar PSF, however when an alternate stellar PSF was used and the 
background was over-subtracted to its uncertainty limit, it was not preferred over the PSF- 
only fit with sufficient confidence, and so is considered unresolved. 

0820+225: The host galaxy of this source is the best resolved of the sample thanks 
to the relatively faint nucleus. Although the host has a strong bulge component, a model 
including both bulge and disk components provides a significantly better fit than the bulge- 
only model if effective radius of the bulge is restricted to < 20 kpc. A compact companion 
galaxy is also resolved at ~0"7 (5.5 kpc) to the west of the nucleus. 

1144S79: This host galaxy is resolved, although not sufficiently well to distinguish 
between a bulge- or disk-dominated morphology. A very faint, arc-like wisp is resolved at 
~2'.'25 (18 kpc) to the northwest of the nucleus. 
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1249+174: This host galaxy is definitely resolved; however, there was insufficient infor- 
mation to constrain its morphology. A bright, elongated structure is resolved extending from 
~0'.'8 to ~1'.'8 (~5 to 13 kpc) to the northwest of the nucleus, while two compact objects are 
visible at ~ 10" from the BL Lac. 

1308+326: The host galaxy of this source is well fitted by an exponential disk profile. 
However, it cannot be confidently distinguished from the PSF if the background is over- 
subtracted to its uncertainty limit, and so it is considered unresolved. Two bright compact 
companion objects are within 6" (48 kpc) of the BL Lac, while a number of fainter objects 
are within half that distance. 

1422+580: This host galaxy is well-resolved, and a de Vaucouleurs profile provides a 
significantly better fit than an exponential disk profile. A host galaxy model including both 
bulge and disk components provides a marginally better fit. The host galaxy appears smooth 
and is clearly elongated, with an axial ratio of 0.8. Several compact companion objects are 
in the near vicinity of the BL Lac, some faint ones within ~ 1'.'5 (13 kpc). 

1533+535: Another well-resolved, smooth-looking host galaxy, although in this case the 
morphology could not be constrained. This may be due to the presence of a disk component, 
as the best-fit de Vaucouleurs effective radius is unphysically large. A de Vaucouleurs + ex- 
ponential profile gives a marginally better fit than the de Vaucouleurs-only model. 

2240-260: The host galaxy of this source is unresolved. A few compact objects lie to 
the northeast, ~ 6" (45 kpc) from the BL Lac. 

4. Host Galaxy Evolution 

To study the evolution of BL Lac host galaxies we combine the STIS observations 
described above (hereafter referred to as the high-z sample) with the lower-redshift BL Lacs 
observed in the WFPC2 snapshot survey (those with z < 0.6; hereafter the low-z sample). 
The results of the host galaxy analysis for these 63 low-z BL Lacs are listed in Table 4, and 
the details of the analysis described in Scarpa et al. (2000). 

The STIS F28x50LP filter samples the rest-frame U to V bands at redshifts z > 0.6, 
and so is particularly sensitive to star formation. This fact has enabled us to detect BL Lac 
host galaxies at greater distances than previously achieved. However, the star formation his- 
tories of the high-z BL Lacs are unknown, and thus the filter conversions and K corrections 
are highly uncertain. Any measurement of luminosity evolution in absolute magnitude space 
would be subject to large errors. However, from NICMOS imaging of 10 low redshift BL Lacs 
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(Scarpa et al. 2000a), we know that these host galaxies have colours consistent with old stel- 
lar populations. At the same time, conversion from F702W to F28x50LP magnitudes at 
lower redshifts is far less sensitive to star formation history than conversion from F28x50LP 
to F702W at higher redshifts. To minimise correction uncertainties we perform the analysis 
of host galaxy luminosity evolution using the measured F28x50LP apparent magnitudes of 
the high-redshift sample, corrected only for galactic extinction, and the derived F28x50LP 
apparent magnitudes of the low-redshift sources, converted from the measured F702W ap- 
parent magnitudes. The filter conversions used for the low-z sample are not fixed, but rather 
varied within the fitting process (described in §4.1). 

The host galaxy brightnesses of the combined sample are then compared to models 
of galaxy luminosity evolution. These models are derived using the GISSEL evolutionary 
synthesis models of Bruzual & Chariot (1993) with a Salpeter IMF (Salpeter 1955), with 
the Padova 2000 library of isochrones (Girardi et al. 2000), and the BaSeL standard stellar 
library (Lejeune, Cuisinier & Buser 1997). 

Figure 5 shows the STIS F28x50LP Hubble diagram for the host galaxies, with the low-z 
sample filter conversions performed assuming a solar metallicity burst population calibrated 
to match the typical early-type colour of B — V = 0.96 (Fukugita, Shimasaku & Ichikawa 
1995). Also shown is the line for a non-evolving population with the same early- type spec- 
trum (dashed line), and a non-evolving late-type spectrum (solar metallicity, age = 8 Gyrs, 
star formation rate decay exponent r = —3) matching the Hubble type Sab colours of 
B — V = 0.57 (Fukugita, Shimasaku & Ichikawa 1995; dotted line). The high redshift host 
galaxies appear to be more luminous than predicted by the non-evolving models. However, 
the large error bars and the presence of non-detections in both the high-z and low-z samples 
necessitate rigorous testing. 



4.1. x 2 Fitting with Censored Data 

The presence of host galaxy non-detections in some of the higher redshift objects means 
that the apparent luminosity evolution may just be a selection effect. Great care must be 
taken in dealing with the non-detections, and for this we employ survival analysis techniques. 

The technique developed to apply survival analysis to the fitting of evolution models 
is derived from the EM algorithm (Dempster, Laird & Rubin 1977). This algorithm uses 
maximum likelihood to estimate the regression coefficients of a 2-dimensional data set, and 
assumes that the censored data points have a normal distribution about the regression line. 
This regression line is determined iteratively, with the estimates for the weighting of the 
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censored data points updated at each iteration. The assumption of a normal distribution 
makes the EM method less general than other, non-parametric methods, but given that the 
host galaxy luminosities have a roughly-normal distribution at low redshifts, this assumption 
seems reasonable. 

As we are fitting the apparent magnitudes rather than the absolute magnitudes, we can't 
perform a simple regression analysis. Instead of iterating over regression line estimates, we 
iterate over evolution models. If we assume that each galaxy consists of a single, co-evolving 
population, then each model is defined by two parameters: magnitude zero point and either 
population age or formation redshift (zf orm ) for the non-evolving and passively evolving 
models respectively. 

Each iteration in population age or formation redshift implies a set of filter corrections 
for the low-z sample. Although there are only small differences in the filter correction between 
different galaxy models at low redshift, we wish to calculate x 2 given the measured F702W 
magnitudes, not the inferred F28x50LP magnitudes, for each evolution model. With this 
in mind, we apply a new, self-consistent spectral correction for each new population age 
or formation redshift, rather than assume the same correction for all iterations. For the 
moment, we assume solar metallicity in all models. 

4-1.1. Single-Population Passive and Non-Evolving Models 

The passively evolving stellar population yielding the best fit to the F28x50LP apparent 
magnitudes of our sample has a formation redshift of Zf orm = 1.8, with a reduced x 2 °f 1-37. 
The la uncertainties are +0.7/ — 0.3 and the 3a uncertainties are +1.0/ — 0.5. A better fit 
by a model with high formation redshift can be ruled out with high certainty — Zf orm = 3.0 
yields a reduced x 2 °f 1-98, while Zf orm = 5 gives a reduced x 2 °f 2.35. An F-test (72 data 
points and 2 degrees of freedom) indicate a worse fit than the Zf orm = 1.8 model with 95% 
and 99% significance respectively. 

The data are also reasonably well fit by an extremely young non-evolving population. 
The best fit is for a population of age 0.02 Gyr, with reduced x 2 = 1-51 — a worse fit than 
the passive model, although without statistical significance. A population as young as 1 Gyr 
provides a poor fit to the data, with reduced x 2 = 2.13: worse than the best passive model 
at the 97.5% significance level, and worse than the 0.02 Gyr model at the 95% significance 
level. A population of age 2.5 Gyrs gives a reduced x 2 of 2.38; worse than the passive model 
at the 99% level. Similar results are obtained for models of non-evolving late-type galaxy 
spectra with the same B — V colours as the burst models. 
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The young populations needed to fit the data do not agree with the colours observed 
in BL Lac host galaxies at low redshift. From Scarpa et al. (2000a), we know that these 
galaxies have R — H = 2.3 ± 0.3, consistent with an old stellar population. At the blue-most 
limit, these colours imply age ~ 5 Gyrs for a burst population. Fitting the normalization of 
a non-evolving model with age = 5 Gyrs results in a reduced x 2 °f 3.05. An F-test indicates 
that the passive evolution model gives a better fit than the 5 Gyr non-evolving model at the 
99.9% significance level. 

Figure 6 shows the Hubble diagrams with the best-fit passive and non-evolving mod- 
els (top and middle). In each case, the low-z host galaxies are converted from F702W to 
F28x50LP magnitudes assuming the corrections implied by the evolution model plotted with 
a solid line. Visual inspection of the Hubble diagrams shows that the non-evolving model 
predicts galaxies about one magnitude too faint at z > 0.6. Even much younger non-evolving 
models (1 Gyrs; middle, dashed line, fail to reproduce the high redshift colours, although 
an unphysically young model (0.02 Gyrs; middle, dotted line), does fit. Passive evolution 
models with high formation redshifts similarly underestimate the brightness of the z > 0.6 
host galaxies, as can be seen by the Zf orm = 5 model (top, dashed line). 



4-1.2. Dual- Population Evolution Model 

While the above results indicate that the luminosity evolution is better fit by a model 
that includes a rapidly evolving stellar population, it may not be the entire population 
that is evolving. If an old stellar population is already in place, then a fresh burst of star 
formation involving only a fraction of the galaxy's mass may also reproduce the observed 
evolution. We fit the data with a model that comprised a passively evolving stellar population 
superimposed with a second population. This second population remains at constant (low) 
age down to a given redshift, after which it also evolves passively, to simulate a normally 
evolving galaxy with some degree of ongoing star formation, which dies out at some point. 
The model now includes five parameters: magnitude offset, formation redshift of the primary 
population (zf orm ), the point at which the secondary population begins to evolve passively 
(z*), percentage mass of the secondary population relative to the primary, m*, and the initial 
age of the secondary population, age*. 

A wide range of parameters provided fits of similar quality to the best-fit one-component 
model. For a primary population with a high formation redshift (zf orm > 4), a secondary 
population of at least 1% of the mass of the primary (depending on z* and age*) produces 
reasonable fits. For a small secondary component of 1% to 3%, z* can be as low as ~ 0.5 while 
still agreeing with the red colours observed in low-redshift BL Lac host galaxies. The best 
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fit was found for m* of 2%, z* = 0.6, age* = 1 Gyr, and Zf orm = 5, with reduced x 2 — 1-29. 

The bottom panel of figure 6 shows the best-fit dual-population evolution model (solid 
line). The dashed line in this plot shows the evolution of the Zf orm = 5 base population. 

Although we can't distinguish between single- and dual-population models, we can say 
with high confidence that the observed luminosity evolution is significantly better represented 
by models which include star formation activity in the range 0.5 < z < 2.5, which diminishes 
rapidly or ceases altogether below z ~ 0.5, than it is by models in which the galaxies are 
undergoing steady star formation to low redshift, or by models in which all star formation 
activity had ceased by z ~ 2.5. 

4.2. Effects of Metallicity and Dust 

Given our knowledge of the red colours of low redshift BL Lac host galaxies, these 
data show conclusively that higher redshift BL Lac hosts are both significantly bluer and 
significantly brighter than their low redshift counterparts. This change is well modelled as 
evolution of the stellar population, however evolution in metallicity or the presence of dust 
may also play a part. 

4-2.1. Evolution in Metallicity 

The evolution results given so far are for populations with constant, solar metallicities. 
If the galaxies at higher redshift have much lower metallicities, the data can be fitted with 
older stellar populations. We simulate non-solar metallicities using the models of Westera et 
al. (2002). 

Taking the host galaxies with z > 0.6 to have extremely low metallicity — 0.02 solar - 
and those with z < 0.6 to have solar metallicity, the best-fit passively evolving population 
has Zf orm = 2.3 (laerrors), with reduced \ 2 — 1-28, better at 97% confidence than a 
Zform — 5 model (reduced \ 2 = 2.07). The same metallicity evolution allows non-evolving 
populations to fit better. A population with age = 1 Gyr gives reduced \ 2 = 1-99, worse 
than the best-fit passive model with only 95% confidence (instead of 97.5% for a constant, 
solar metallicity). However this population still does not satisfy the low- redshift R — H 
colour constraint. 

Of course, a sudden transition from 0.02 solar to 1 solar metallicity at z — 0.6 is not 
realistic; however, this scenario provides a limit on the extent to which metallicity evolution 
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can affect the results. 



4-2.2. Effects of Dust Reddening on the Fits 

The presence of large amounts of dust in the host galaxies could potentially result in 
extreme reddening — enough to give a young population the R — H colours observed in low 
redshift BL Lacs. However, there is a limit to how much dust can be added, as the preferential 
extinction caused by dust at shorter wavelengths makes the high-redshift population fainter, 
worsening the fit. 

We applied models of dust extinction from Ferrara et al. (1999) to our spectral evolution 
models to determine whether a dust-reddened, non-evolving population could produce a 
reasonable fit to the data, with the constraint that the low- redshift (z < 0.4) host galaxies 
have R — H = 2.0 (the bluest possible based on the NICMOS observations), after dust 
reddening. The extinction model which achieved this most efficiently was one based on the 
Milky Way extinction curve, assuming an elliptical galaxy geometry with a Jaffe density 
profile, and a central V-band optical depth of 7y(0) = 2. With this model, a non-evolving, 
1 Gyr-o\d population has R — H = 2 at low redshift, and fits our Hubble diagram with 
reduced x 2 — 2-21. This is worse than the best-fit passive model at the 98% level by the 
F-test. To substantially improve this fit (worse only at the 95% level), a 0.6 Gyr population 
is needed, which, with this dust model, has R — H = 1.7 — certainly far too blue for the 
low-redshift host galaxies. 

Thus, no (non-evolving) dust model can make a non-evolving stellar population fit the 
data. The level of dust reddening needed even for these relatively poor fits, ry(0) = 2, is 
at the upper limit of dust optical depths estimated in ellipticals by Wise & Silva (1996) or 
face-on spirals by Kuchinski et al. (1998), and is substantially higher that the limit predicted 
for ellipticals by Goudfrooij & de Jong (1995) of ry(0) < 0.7, who find more typical values 
of 7y(0) -0.2-0.3. 

4-2.3. Evolution in Dust Content 

It is possible to reproduce the observed luminosity evolution for a non- or weakly- 
evolving population if dust content is allowed to increase rapidly between z ~ 0.8 and ~ 0.4. 
Large changes at lower redshifts fail to fit the weak evolution seen at low-redshift, and 
extending the dust evolution to higher redshifts means that there is too much attenuation 
to reproduce the bright host galaxies at z > 0.6. 
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Assuming that the galaxies are relatively dust-free at z > 0.8, and that optical depth 
in dust increases linearly to z — 0.4, the minimum final optical depth required to produce a 
reasonable fit to the data and to satisfy the low-z R — H = 2 constraint is 7y(0) = 8. This 
allows a non-evolving model as young as 0.8 Gyrs to fit the data with reduced x 2 = 1-54, 
within la - of the best-fit passive evolution model. 

Of course, Ty(0) = 8 is extremely high — over an order of magnitude higher than the 
typical values determined for local ellipticals, and several times higher than those determined 
for star-forming galaxies (e.g. Calzetti 2001). If we assume a local optical depth of Ty(0) = 2, 
the upper limit predicted for ellipticals, then the best-fit non-evolving evolution model gives 
reduced x 2 — 2.01 — still worse than the best passive model at the ~95% level. 

Approximately same level of dust evolution is required to make a passively evolving 
population with high formation redshift fit the data. A population with Zf orm = 5 requires 
a local optical depth of 7y(0) = 10 for reduced x 2 — 1-52, while the same population with 
7y(0) = 2 at z = gives reduced x 2 — 2-09 — worse than the -2/ orm = 1.8 model at 97% 
confidence. 

A concurrent evolution in metallicity will allow these models to fit the data with a smaller 
increase in dust content. Using the metallicity evolution described above, the minimum 
local optical depth needed for a 2 Gyr-old non-evolving or a Zf orm = 5 passively evolving 
population to fit the data (within la of the best passive model) is t v (0) = 5 (reduced x 2 ~ 
1.5). This is still an extreme increase in dust content. If 7y(0) = 2 locally, then these models 
are worse than the best passive model at the ~95% confidence level. 

Figure 7 shows the Hubble diagram with both the non-evolving and passive evolution 
models, dust and metallicity evolution included. 

4.3. Absolute Luminosity Evolution 

The lack of extended emission line light or a big blue bump has made it possible to use 
the extremely sensitive F28x50LP filter in this study. However, to compare our evolution 
results directly to other studies, it is necessary to determine this evolution in absolute near- 
infrared magnitudes. This result depends strongly on the K corrections used. To find the full 
range of possible magnitude-redshift gradients, we perform regression analysis on absolute 
K-band magnitudes using the range of K corrections implied by the full range of possible 
evolutionary models from the above analysis. 

We use the EM algorithm to perform regression analysis while accounting for the cen- 
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sored data points, via the IRAF task EMMETHOD. This task by itself does not take error 
bars into account. As the uncertainties in host galaxy magnitudes increase with redshift, 
they are an important consideration. To account for them, we use a Monte Carlo approach. 
One hundred simulated data sets were generated, randomizing each detected data point over 
a composite normal distribution, with the widths of the upper and lower halves of the dis- 
tribution defined by the upper and lower error values for that data point. The regression 
coefficients and their uncertainties were determined for each of these 100 simulated data set 
with EM method, and the final coefficients taken to be the mean values. The uncertainty 
values were taken to be the mean of the individual uncertainties in the regression coeffi- 
cients, added in quadrature with the standard deviation of the coefficients across the 100 
simulations. 

The steepest slope occurs when K corrections are derived from the evolutionary model 
with the oldest possible population at high redshift (at the 3a level, assuming no dust or 
metallicity evolution) - - the passive evolution model with Zf orm = 2.8. The shallowest 
slope is in the case of the youngest possible high-z population (also at the 3a level) — a 
passive model with Zf orm = 1.3. Unsurprisingly, the best-fit passive evolution model, with 
z form — 1-8, gives a gradient about halfway between these. The results of the regression 
analysis using these three sets of K corrections are shown in Table 5. 

A similar range of gradients are obtained if we assume a secondary burst component 
added to an older population, with a lower limit of 0.59 mag/z for a passive, old population 
plus a 0.1 Gyr burst involving 3% of the galaxy's mass which itself passively evolves below 
z = 0.5. For a constant-age, old population, or a population with high formation redshift 
we obtain ~ 1.7 mag/z. To obtain high redshift K-band magnitudes consistent with no 
K-band luminosity evolution, the stellar populations at z > 0.6 must have ages less than 
0.2 Gyrs. Such star formation histories can only be made to fit the observed Hubble diagram 
and satisfy the low-redshift colour constraints with extreme levels of dust evolution (from 
7y(0) ~ to 7y(0) ~ 5 between z = 0.8 and 0.4). 

Using the best fit to the Hubble diagram to derive K corrections, and adding in quadra- 
ture the error in the regression analysis to the probable range of slopes from different K 
corrections, we determine the overall K-band luminosity evolution of BL Lac host galax- 
ies to be 1.0 ± 0.45 mag/z. Figure 8 shows absolute K-band magnitude versus redshift, K 
corrected assuming a passively evolving population with Zf orm = 1.8. Also plotted are the 
best-fit linear regression lines representing the magnitude-redshift gradients derived assuming 
the three different passive evolution models. 
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5. Discussion 

We detect luminosity evolution in BL Lac host galaxies, consistent with the primary 
stellar population evolving strongly from 1.5 < z < 2.5, or with a smaller burst of star 
formation evolving from z > 0.5, and inconsistent with a non-evolving population or a 
population evolving passively from high redshift. Interestingly, this formation epoch for 
the primary population corresponds roughly to the epoch of maximum quasar activity (e.g. 
Maloney & Petrosian 1999; Shaver et al. 1999). This is consistent with a scenario in which 
black hole feeding and growth — i.e., the epoch of maximum AGN activity — coincides with 
the formation of the surrounding bulge (Kormendy & Gebhardt 2001). 

Radio galaxies (Djorgovski, Spinrad & Marr 1985; Spinrad 1986; Spinrad & Djorgovski 
1987; Lilly & Longair 1982; Lilly, Longair & Allington-Smith 1985; Ar agon- Salamanca et 
al. 1993; McLure & Dunlop 2000) and the hosts of radio-loud quasars (Lehnert et al. 1999; 
Kukula et al. 2001; Hutchings 2001) exhibit similar evolution, probably dimming by of order 
1 magnitude in K between z = 1 and 0. Although this evolution is not as well constrained, 
it is consistent with that seen in our sample. This is in agreement with the unified model 
of AGN in which radio galaxies, radio quasars and blazars are expected to have the same 
host galaxy properties. In contrast, BCGs (Aragon-Salamanca et al. 1993) other quiescent 
early-types (Lilly et al. 1995; Stanford et al. 1998; Bell et al. 2004) and radio-quiet quasar 
hosts (Kukula et al. 2001; Rix et al. 2001; Ridgway et al. 2001) exhibit flat, or even negative 
luminosity evolution, with local galaxies appearing as bright or brighter than those at z > 1. 
Thus, the hosting of a powerful radio jet seems to be linked to a galaxy's luminosity evolution. 

The luminosity evolution of BL Lac host galaxies is unambiguous, not affected by the 
biases inherent in other radio-loud classes (see §1). We can therefore ascribe the evolution 
observed in BL Lac hosts to a real change in galaxy luminosities, which in turn grants us 
more confidence that the evolution observed in other radio-loud hosts is real. 

Furthermore, BL Lacs allow us significantly more confidence that the observed evolution 
is a luminosity evolution in individual host galaxies. There is mounting evidence to suggest 
that BL Lac objects are long-lived, possibly with lifespans on the same timescales as galaxy 
evolution. The most compelling is the flat, possibly negative number density evolution of 
many BL Lacs (Morris et al. 1991; Giommi, Menna & Padovani 1999; Bade et al. 1998; 
Perlman et al. 1996). There are as many, if not more BL Lac objects now than during the 
epoch of peak quasar activity. The low accretion rates in BL Lacs (less than 1% Eddington; 
O'Dowd, Urry & Scarpa 2002; Wang, Staubert & Ho 2002) also suggest that BL Lacs may 
burn stably for extended periods. If the jets of BL Lacs obtain most of their power from the 
rotational energy of the black hole (e.g. Blandford & Zjanek 1977, Meier 2002), they are 
expected to have lifespans of the order of gigayears (Cavaliere & Malquori 1999). Long-lived 



- 21 - 



emission is also supported by their dynamically evolved state at low redshift (e.g., Urry et 
al. 2000) compared to the high frequency of interaction remnants and close companions at 
higher redshifts (see §3.3). Thus, with BL Lacs we may be observing luminosity evolution 
intrinsic to the host galaxies themselves, rather than evolution in the processes which select 
galaxies as hosts of radio-loud AGN. 

Our results are consistent with and confirm the results of the ground-based study of 
Heidt et al. (2004), who find evidence of luminosity evolution in BL Lac host galaxies, 
although do not constrain this evolution. Of the sources which overlap the two samples - 
0820+225, 1249+174, 1422+580 and 2240-260 — Heidt et al. detect the host of 0820+225. 
Our host galaxy magnitudes for this source agree within error, although our effective radii 
do not. We believed that this stems from the uncertainty and width of the ground-based 
PSF (of order that of the host itself in some cases). 

5.1. Possible Causes of Luminosity Evolution 

Both quiescent, early-type galaxies and radio-quiet AGN host galaxies are gaining lu- 
minous mass up to the present, in line with the predictions of hierarchical clustering models 
(Kauffmann & Haehnelt 2000), therefore the luminosity of these galaxies does not change 
with time. This does not appear to be the case with the host galaxies of radio-loud AGN, 
whose dimming over cosmic time is apparently not offset by the accretion of significantly 
amounts of new material. A number of effects may contribute to this. 

Evolution linked to emission processes: Low redshift radio-loud AGN used in evolution 
studies are dominated by sources with low accretion onto their SMBHs, such as the FR I 
and low-power FR II radio galaxies in the 3C sample (in which luminosity evolution has 
been best studied). These are also the proposed parent population of BL Lacs (Urry & 
Padovani 1995). If a large fraction of these AGN are long-lived, jet-dominated sources, then 
we select against active recent merger histories when we observe them locally. Continued 
merger events following an AGN's formation epoch will refuel the nucleus after depletion of 
the initial supply of gas and possibly disrupt stable, jet-dominated emission (Cavaliere & 
Padovani 1989), as well as add luminous mass to the galaxy, offsetting the dimming due to 
an aging primary stellar population. 

Evolution Linked to Black Hole Properties: While the production of a radio jet is 
likely to be linked to the properties of the central SMBH, it is uncertain exactly which black 
hole properties are important in determining radio loudness. There is some evidence that 
high black hole mass plays a part (Franceschini, Vercellone & Fabian 1998; Lacy et al. 2001; 
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Nagar et al. 2002; Jarvis k McLure 2002; O'Dowd, Urry & Scarpa 2002; Dunlop et al. 
2003), although this may just be a selection effect (Woo & Urry 2002; Oshlack, Webster 
& Whiting 2002; Ho 2002). Black hole spin seems the more likely candidate as the critical 
determinant (Meier 2002; O'Dowd et al. 2004). On the other hand, high host galaxy mass 
is expected to dampen radio emission (Bicknell 1995). Galaxies that have extreme SMBH 
properties relative to their mass should favor the production of radio jets. Such SMBHs 
may be expected to preferentially occur in those galaxies whose formation was completed by 
higher redshifts. 

In the case of SMBH mass, hierarchical models (Kauffmann & Haehnelt 2000) predict 
that earlier-formed galaxies will have higher SMBH-to-host galaxy mass ratios, due to the 
larger amounts of cold gas available for accretion at earlier epochs. A correlation has even 
been observed between the fraction of a galaxy's mass in the black hole and the age of its 
stellar population (Merrifield et al. 2000). 

In the case of spin, a black hole of a given mass that is formed from a small number of 
major merger events is likely to have a higher spin than one of the same mass formed from 
many less massive mergers, as the angular momenta of these lesser interactions will tend to 
cancel out unless they are correlated (Cavaliere & Padovani 1989). Locally, this means black 
holes of given mass which formed during the epoch of major mergers, at z ~1 - 2, will have 
higher spin than whose whose formation proceeded to the present. 

According to the scenarios outlined above, when we select radio sources at low redshift, 
we are dominated by low-accretion AGN with extreme SMBH properties that are prefer- 
entially produced at higher redshifts, and that should have had relatively quiescent merger 
histories since the epoch of formation of their central engines. We have more credible evi- 
dence that radio-loud AGN do exhibit just such merger histories, supporting the idea that 
the properties of radio- loud AGN SMBHs are typically defined at high redshift. 

5.2. The Merger History of Radio Hosts 

In §3.3 it was found that many of the BL Lacs in the high-z sample have either close 
companions (62.5%) or signs of recent interactions (50%). Fewer of the broader sample of 
Urry et al. (2000) had close companions (47%) and only a handful (<10%) showed signs of 
interactions. This supports the idea that the high-z BL Lacs are closer to their formation 
epoch, and that the observed luminosity evolution is linked to the merger history of these 
host galaxies. 

It is impossible to perform a such a comparison over a similar redshift range for other 
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radio-loud AGN, as flux limits induce a power-redshift correlation. Locally, however, FR I 
radio galaxies, which are thought to be the primary parent population of BL Lac objects, 
do show fewer signs of recent interactions than the more powerful FR lis (Heckman et al. 
1986; Smith & Heckman 1989), with the latter hosted by bluer galaxies and showing a much 
higher incidence of distorted morphologies. 

These observations, along with the scenarios outlined in §5.1, point to a scheme in 
which radio-loud AGN are formed in major mergers, typically at high redshift, while those 
observed locally have evolved from a significantly earlier formation epoch relatively free of 
further interactions. 

Given this picture, it is expected that cluster environment will play a critical role in 
governing when and where radio-loud AGN form and survive. Formation will be dictated by 
a balance between a high density environment (clusters for merger frequency) and low cluster 
velocity dispersion (high redshift clusters or local field, to permit full mergers; Makino & 
Ebisuzaki 1996). Survival as jet-dominated sources to low redshift will favour lower density 
environments, as we require fewer ongoing merger events. BL Lacs do prefer relatively poor 
cluster environments at low redshift (z < 0.4), of Abel richness class ~0 to 1 (Falomo, 
Pesce & Treves 1993; Pesce, Falomo & Treves 1994; Fried, Stickel & Kiihr 1993; Wurtz et 
al. 1997; Smith, O'Dea & Baum 1995). At higher redshifts, their environments are not 
well known, although at z > 0.6 at least some are found in rich clusters (Fried, Stickel 
& Kiihr 1993; Wurtz et al. 1997). This is consistent with the picture in which BL Lacs 
are preferentially produced in cluster environments at high redshift (when cluster velocity 
dispersions were lower), but survive only in the sparcer regions of these clusters. Reports 
on the environmental evolution of other classes of radio-loud AGN are conflicting, as are the 
environmental differences between radio-loud and radio-quiet AGN. 



6. Conclusions 

Nine BL Lac objects with z > 0.6 were imaged with HST, using the STIS CCD and 
the F28x50LP filter. The host galaxies of five of these were resolved with high confidence. 
For all but 0820+225, these are the first detections of host galaxies in these sources, and 
include the two highest redshift BL Lac host galaxies detected to date. These galaxies are 
bright, with an average absolute magnitude of <M R >= —24.3 mag assuming an early-type 
spectrum and —23.7 mag assuming a late-type spectrum. 

The morphologies of two of the galaxies were well resolved. Both 1422+580 and 0820+225 
are better fit by a de Vaucouleurs profile than by an exponential disk profile, indicating bulge 
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dominated hosts, as is typical for radio-loud AGN. However, the host galaxy of 0820+225 
was significantly better fit by a composite model including both bulge and disk components 
than by a bulge component alone. Many, including some of those technically considered 
non-detections, show evidence of having extended, disky components. A large fraction of 
the sample show evidence of close companions and distorted morphologies — far more than 
are observed at lower redshift. 

Combining the derived host galaxy magnitudes with the low redshift host galaxies from 
Urry et al. (2000), we measure luminosity evolution in a sample of 72 HST-imaged BL Lac 
host galaxies spanning the range < z < 1. We find that the Hubble diagram is well fit by 
a model of a stellar population evolving passively from z — 1.8 1 ' 3 . Alternatively, an older 
population with a small fraction (<3%) undergoing active star formation to z ~ 0.5 also 
fits the data well. A single population undergoing passive evolution from a high formation 
redshift (zf orm > 5) can be ruled out with high confidence. By constraining the colours of the 
low redshift host galaxies based on NICMOS observations, models including no luminosity 
evolution can also be ruled out with high confidence. Non-evolving and high-zj orm passive 
evolution models are still ruled out when reasonable models of dust and metallicity evolution 
are included. 

By K-correcting the host galaxy magnitudes based on the range of acceptable fits to 
the Hubble diagram, we determined an average K-band luminosity evolution of 1.0 ± 0.45 
magnitudes per unit redshift. 

The fact that flux selection of BL Lacs is dominated by beaming angle rather than 
intrinsic brightness, together with the low levels of extended emission from ionized gas and 
scattered nucleur light in these sources, mean that they are relatively free of the biases af- 
fecting quasar studies. This increases our confidence that the observed luminosity evolution 
is real, and in turn in the luminosity evolution observed in other classes of radio-loud AGN. 
We also have reason to believe that the observed evolution is intrinsic to the galaxies them- 
selves, as evidence suggests that BL Lacs may have lifespans comparable to the timescale 
for galaxy evolution. 

The luminosity evolution of BL Lac host galaxies is steeper than that observed in radio- 
quiet host galaxies, brightest cluster galaxies, and quiescent early-type galaxies, and is con- 
sistent with that observed in other radio-loud hosts. The absence of luminosity evolution in 
radio-quiet and quiescent galaxies suggests that dimming due to an aging population is offset 
by the accretion of new material, as is predicted by hierarchical clustering models. This ap- 
pears not to be the case in the host galaxies of radio-loud sources, which evolve strongly from 
some redshift greater than z ~ 0.5, but also must have experienced significant star formation 
activity in the range 0.5 < z < 2.5. This is further supported by the fact that BL Lacs with 
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z > 0.6 are much more frequently in interacting systems, and have more close companions 
than those at low redshift, which appear dynamically evolved. These observations lead to 
the following conclusions: 

First, it seems that the physical conditions that allow a galaxy to generate a radio jet 
are preferentially produced in a particular cosmic epoch, albeit a broad one: 0.5 < z < 2.5. 
This may be linked to the amount of gas available for fueling, which is also believed to define 
the quasar epoch. However, to explain the radio-loud — radio-quiet dichotomy, it must also 
be linked to the production of the right SMBH conditions. 

Second, radio-loud AGN hosts at low redshift accrete less new material after the epoch 
in which they seemingly form their SMBHs than do other early-type galaxies. A number of 
effects may result in the selection of inactive recent merger histories when we select radio- 
loud AGN at low redshift; however, all point to a scenario in which many radio-loud AGN 
observed at low redshift have survived from a much earlier formation epoch. 

The evolution observed in the stellar populations of radio-loud AGN suggests the black 
holes which power these radio sources are born in gas-rich interactions at z > 0.5. At low 
redshift, where samples are dominated by low-accretion sources, these observations suggest 
that we are observing the last embers of the radio quasar epoch — old AGN with massive, 
rapidly-rotating black holes which have experienced little merger activity since the formation 
of their central engines. 

Support for this work was provided by NASA through grant numbers HST-GO-09121.01- 
A and HST-GO-09223.01-A from the Space Telescope Science Institute, which is operated 
by AURA, Inc., under NASA contract NAS5-26555. 
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Table 1. 



Name 


z 


Type a 


m^(nucl.) b 


M R (nucl.) c 


m^(host) d 


Exp. Time (s) 


Obs. Date) 


0138-097 


0.733 


L 


17.7 


-25.57 


>19.57 


2554.0 


24/07/2001 


0235+164 


0.94 


L 


18.90 


-25.03 


>19.4 


5429.0 


18/07/2001 


0820+225 


0.951 


L 


20.0 


-23.97 


>21.25 


5473.0 


17/09/2001 


1144-379 


1.048 


L 


18.0 


-26.22 


>22.72 


5615.0 


04/08/2001 


1249+174 


0.644 


H 


18.5 


-24.43 


>21.3 


2552.0 


21/07/2001 


1308+326 


0.997 


L 


18.1 


-25.99 


>20.13 


5551.0 


22/07/2001 


1422+580 


0.683 


H 


19.07 


-24.01 


>21.7 


4644.0 


26/08/2001 


1533+535 


0.89 


H 


18.60 


-25.19 


>19.45 


4861.0 


14/09/2001 


2240-260 


0.774 


L 


17.5 


-25.91 


>21.45 


2578.0 


27/07/2001 



a SED type: H=HBL (F lkeV /F 5GHz > 5.5) and L=LBL (F lkeV /F 5GHz < 5.5) 
b Apparent R-band magnitude of source, not extinction or K corrected. 

c Absolute R-band magnitude of source, extinction corrected and K corrected assuming a — — 1, 
H = 70 km s _1 Mpc _1 , Vt M = 0.3, Q x = 0.7. 

d imit on apparent R- magnitude of the host galaxy from Urry et al. (2000). 



Table 2. Results of two-component host galaxy fits 



Name 


Modcl a 


mj^p nucleus^ 


M R nucl. c 


m^p host b 


Mr host c 




(arcsec 


d 


re (kpc) d 


b/a e 


x 2f 




non-detect. 


16.63 


— 0.01/+0.01 


-26.78 


>19.19 




> -25.31 












0.971 


0138-097 


(de Vauc.) 


M 6 74"t 




(-26.67) 


(20.01) 




(-24.33) 


(4.2) 






(30.6) 




0.959 




(cxpon. ) 


(16.69) 




(-26.72) 


(20.36) 




(-23.52) 


(1.1) 






(8.0) 




0.962 




non-detect. 


17.93 


— 0.01/+0.02 


-26.14 


>20.61 




> -25.04 












1.055 


0235+164 


(de Vauc.) 


(18.09) 




(-25.98) 


(20.80) 




(-24.69) 


(2.4) 






(18.9) 




1.022 




(cxpon. ) 


(17.97) 




(-26.10) 


(21.16) 




(-23.69) 


(1.9) 






(15.0) 




1.018 




dc Vauc. 


20.38 


— 0.13/+0.11 


-23.81 


21.16 


—0.25 /+0.40 


-24.64 


5.7 


f 14.2/- 


-5.4 


45.2 


1.0 


0.847 


0820+225 


expon. 


20.24 


-0.12/+0.16 


-23.95 


21.61 


-0.80/+0.84 


-23.54 


1.0 


+0.3/- 


-0.4 


7.9 


1.0 


0.886 




PSF only 


20.28 


-0.05/+0.04 


-23.91 


















0.928 




de Vauc. 


18.00 


-0.12/+0.06 


-26.74 


21.18 


-0.46/+0.38 


-25.47 


2.5 


+4.0/- 


-0.9 


20.2 


1.0 


1.521 


1144-379 


expon. 


17.99 


-0.09/+0.09 


-26.75 


21.28 


-0.42/+0.20 


— 24.52 


0.9 


+0.5/- 


-0.3 


7.3 


1.0 


1.547 




PSF only 


17.81 


-0.02/+0.03 


-26.93 


















1.564 




dc Vauc. 


18.94 


-0.12/+0.08 


-24.09 


21.22 


-0.25/+0.62 


-22.71 


1.1 


+ 2.3/- 


-0.7 


7.6 


1.0 


0.452 


1249+174 


cxpon. 


18.68 


-0.11/+0.08 


-24.35 


22.23 


-0.24/+0.52 


-21.29 


0.8 


+ 1.7/- 


-0.6 


5.5 


1.0 


0.465 




PSF only 


18.72 


-0.04/+0.05 


-24.31 


















0.488 




non-detect. 


17.66 


-0.04/+0.04 


-26.47 


>21.29 




> -24.59 












0.776 


1308+326 






























(cxpon. ) 


(17.69) 




(-26.44) 


(21.45) 




(-23.62) 


(3.0) 






(24.0) 




0.762 




dc Vauc. 


18.56 


-0.08/+0.06 


-24.59 


20.35 


-0.30/+0.22 


-23.78 


2.6 


+ 1.4/- 


-1.1 


18.4 


0.8 ±0.09 


1.074 


1422+580 


cxpon. 


18.35 


-0.04/+0.06 


-24.80 


20.88 


-0.10/+0.06 


-22.81 


0.9 


+ 1.1/- 


-0.5 


6.4 


0.75 ±0.13 


1.140 




PSF only 


18.16 


-0.03/+0.03 


-24.99 


















1.296 




dc Vauc. 


18.26 


-0.07/+0.06 


-25.59 


20.63 


-0.24/+0.33 


-24.66 


3.4 


+3.0/- 


-1.0 


26.4 


1.0 


0.517 


1533+535 


cxpon. 


18.18 


-0.08/+0.08 


-25.67 


21.13 


-0.50/+0.31 


-23.61 


0.8 


+0.6/- 


-0.2 


6.2 


1.0 


0.541 




PSF only 


17.99 


-0.02/+0.04 


-25.86 


















0.624 


2240-260 


non-detect. 


16.64 


-0.02/+0.02 


-26.84 


>20.58 




> -24.08 












0.598 



co 

NO 



a Model fitted: PSF+de Vaucouleurs, PSF+exponential disk, or PSF only; non-detect indicates fit with PSF alone where host galaxy is not resolved. Bracketed 
values indicate best host galaxy fits, although arc technically non-dctcctions. 

^Apparent magnitude of nucleus/host galaxy in STIS F2Sx50LP band, not extinction corrected or K corrected. 

c Absolute magnitude of nucleus/host galaxy in Cousins R band, corrected for galactic extinction and K corrected assuming a — 1 power law for nuclei, early-type 
spectrum for de Vaucouleurs hosts and late-type spectrum for exponential disk hosts, and Hq — 70 km s — 1 Mpc~ 1 , fijvf — 0.3, €l\ — 0.7. 

■^Characteristic radius, assuming g — for kpc radii. 

°Axial ratio. 

f Reduced x 2 f° r fit, calculated comparing full two-dimensional model to the image. 
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Table 3. Results of three-component host galaxy fits 



Name 


M H (nucl.) a 


mji(nvLcl.) b 


Component 


rri£p(host) d 


M fl (host) e 


r e (") f 


r e (kpc)s 


x 2h 


0820+225 


20.37 


-25.14 


de Vauc. 


22.18 


-24.85 


1.4 


11.0 


0.811 








expon. 


21.63 


-24.75 


2.0 


15.8 




1144-379 


17.98 


-28.12 


de Vauc. 


21.35 


-26.28 


1.1 


8.9 


1.514 








expon. 


21.90 


-24.88 


1.3 


10.5 




1249+174 


18.71 


-25.54 


de Vauc. 


21.91 


-23.37 


0.7 


4.8 


0.446 








expon. 


22.79 


-22.08 


1.9 


13.1 




1422+580 


18.37 


-26.02 


de Vauc. 


21.00 


-24.82 


1.4 


9.9 


1.049 








expon. 


21.33 


-24.05 


2.4 


17.0 




1533+535 


18.25 


-26.90 


de Vauc. 


21.10 


-25.52 


1.8 


14.0 


0.505 








expon. 


21.59 


-24.47 


3.0 


23.3 





a Apparcnt magnitude of nucleus in STIS F28x50LP band, not extinction or K corrected. 

b Absolutc R-band magnitude of nucleus, corrected for galactic extinction and K corrected assuming a = 1, 
H = 70 km s^Mpc- 1 , Q M = 0.3, Q x = 0.7. 

c de Vauc. = de Vaucouleurs (bulge) component; expon. = exponential disk component. 

d Apparent magnitude of host galaxy in F28x50LP band, not extinction or K corrected. 

° Absolute R-band magnitude of host galaxy, extinction corrected and K corrected assuming early- and late-type 
spectra for de Vaucouleurs and exponential components respectively, with Ho = 70 km s _1 Mpc _1 , Dm = 0.3, 
U x = 0.7. 

f Characteristic radius in arcseconds. 

g Characteristic radius in kiloparsecs, assuming go = 0. 

h Reduccd \ 2 f° r fit, calculated comparing full two-dimensional model to the image. 
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Table 4. The low-z sample: BL Lac host galaxies from Urry et al. (2000). 



Name 


z 


Type a 


mn nucleus 13 




host c 


mji nucl. d 


niR host c 




arcsec) 


r e (kpc) s 


0118-272 


0.559 


L 


15.78 


± 


0.10 


>19.09 




-26.78 


> -23.47 




— 




0.0 


0122+090 


0.339 


H 


21.98 


± 


0.25 


18.88 


± 


0.04 


-19.29 


-22.85 


1.05 


± 


0.10 


5.08 


0158+001 


0.229 


H 


18.38 


± 


0.06 


18.27 


± 


0.03 


-21.91 


-22.29 


1.90 


± 


0.10 


6.96 


0229+200 


0.139 


H 


18.58 


± 


0.35 


15.85 


± 


0.01 


-20.52 


-23.41 


3.25 


± 


0.07 


7.97 


0257+342 


0.247 


H 


19.18 


± 


0.30 


17.93 


± 


0.01 


-21.30 


-22.85 


1.75 


± 


0.12 


6.78 


0317+183 


0.190 


H 


18.28 


± 


0.05 


17.59 


± 


0.01 


-21.56 


-22.46 


3.25 


± 


0.10 


10.31 


0331-362 


0.308 


H 


19.03 


± 


0.10 


17.81 


± 


0.02 


-22.00 


-23.62 


3.10 


± 


0.20 


14.07 


0347-121 


0.188 


H 


18.28 


± 


0.15 


17.72 


± 


0.01 


-21.53 


-22.30 


1.25 


± 


0.05 


3.93 


0350-371 


0.165 


H 


18.03 


± 


0.15 


17.08 


± 


0.01 


-21.46 


-22.60 


1.70 


± 


0.07 


4.81 


0414+009 


0.287 


H 


16.08 


± 


0.05 


17.49 


± 


0.02 


-24.77 


-23.71 


4.70 


± 


0.50 


20.31 


0419+194 


0.512 


H 


19.53 


± 


0.17 


21.05 


± 


0.15 


-22.80 


-22.22 


0.40 


± 


0.07 


2.47 


0502+675 


0.314 


H 


17.33 


± 


0.10 


18.86 


± 


0.09 


-23.75 


-22.62 


0.60 


± 


0.07 


2.76 


0506-039 


0.304 


H 


18.73 


± 


0.15 


18.35 


± 


0.01 


-22.27 


-23.03 


0.60 


± 


0.05 


7.19 


0521-365 


0.055 


L 


15.28 


± 


0.10 


14.60 


± 


0.01 


-21.68 


-22.41 


2.80 


± 


0.07 


3.00 


0548-322 


0.069 


H 


16.93 


± 


0.10 


14.62 


± 


0.01 


-20.54 


-22.92 


7.05 


± 


0.15 


9.31 


0607+710 


0.267 


H 


18.23 


± 


0.10 


17.83 


± 


0.02 


-22.44 


-23.17 


2.40 


± 


0.12 


9.85 


0706+591 


0.125 


H 


17.53 


± 


0.07 


15.94 


± 


0.01 


-21.31 


-23.04 


3.05 


± 


0.07 


6.84 


0735+178 


0.424 


L 


16.58 


± 


0.07 


>20.44 




-25.26 


> -22.06 








0.0 


0737+744 


0.315 


H 


17.88 


± 


0.15 


18.01 


± 


0.08 


-23.20 


-23.48 


2.10 


± 


0.45 


9.68 


0806+524 


0.138 


H 


15.98 


± 


0.02 


16.62 


± 


0.01 


-23.09 


-22.61 


1.45 


± 


0.03 


3.53 


0823+033 


0.506 


L 


17.78 


± 


0.11 


>20.18 




-24.52 


> -23.04 








0.0 


0828+493 


0.548 


L 


18.93 


± 


0.12 


20.26 


± 


0.10 


-23.57 


-23.31 


0.65 


± 


0.10 


4.16 


0829+046 


0.180 


L 


15.88 


± 


0.07 


16.94 


± 


0.04 


-23.83 


-22.98 


4.30 


± 


0.75 


13.06 


0851+202 


0.306 


L 


14.99 


± 


0.06 


>18.53 




-26.02 


> -22.88 








0.0 


0927+500 


0.188 


H 


17.48 


± 


0.30 


17.62 


± 


0.05 


-22.33 


-22.40 


2.00 


± 


0.45 


6.29 


0954+658 


0.367 


L 


16.08 


± 


0.06 


>19.60 




-25.39 


> -22.40 








0.0 


0958+210 


0.344 


H 


21.48 


± 


0.40 


18.93 


± 


0.01 


-19.82 


-22.84 


0.82 


± 


0.04 


4.01 


1011+496 


0.2 


H 


15.88 


± 


0.05 


17.30 


± 


0.02 


-24.08 


-22.89 


1.80 


± 


0.12 


5.94 


1028+511 


0.361 


H 


16.48 


± 


0.10 


18.55 


± 


0.08 


-24.95 


-23.40 


1.80 


± 


0.35 


9.08 


1 1 04-1-384 

J_ 1U4 1 *J O *± 


0.031 


H 


13.78 


± 


0.08 


13.29 


± 


0.02 


—21.90 


—22.41 


3.95 


± 


0.05 


2.45 


1133+161 


0.460 


H 


20.28 


± 


0.18 


19.83 


± 


0.04 


-21.77 


-22.99 


1.55 


± 


0.23 


9.05 


1136+704 


0.045 


H 


16.15 


± 


0.04 


14.45 


± 


0.02 


-20.35 


-22.10 


3.10 


± 


0.02 


2.75 


1212+078 


0.136 


H 


16.38 


± 


0.10 


16.02 


± 


0.01 


-22.66 


-23.17 


3.40 


± 


0.10 


8.19 


1215+303 


0.130 


H 


14.55 


± 


0.01 


15.99 


± 


0.01 


-24.38 


-23.08 


8.35 


± 


0.20 


19.36 


1218+304 


0.182 


H 


15.68 


± 


0.10 


17.12 


± 


0.03 


-24.05 


-22.82 


2.60 


± 


0.30 


7.97 


1221+245 


0.218 


H 


16.89 


± 


0.05 


18.63 


± 


0.06 


-23.28 


-21.79 


1.25 


± 


0.25 


4.42 


1229+643 


0.164 


H 


18.03 


± 


0.30 


16.38 


± 


0.01 


-21.45 


-23.29 


2.00 


± 


0.07 


5.63 


1248-296 


0.370 


H 


18.83 


± 


0.08 


18.87 


± 


0.02 


-22.66 


-23.15 


1.10 


± 


0.05 


5.63 


1255+244 


0.141 


H 


17.08 


± 


0.05 


16.72 


± 


0.01 


-22.05 


-22.57 


2.50 


± 


0.05 


6.21 


1402+041 


0.340 


H 


16.38 


± 


0.01 


>19.38 




-24.89 


> -22.36 








0.0 


1407+595 


0.495 


H 


18.84 


± 


0.05 


19.04 


± 


0.05 


-23.40 


-24.08 


1.75 


± 


0.38 


10.63 


1418+546 


0.152 


L 


15.68 


± 


0.06 


16.10 


± 


0.02 


-23.62 


-23.37 


3.65 


± 


0.11 


9.65 


1426+428 


0.129 


H 


17.38 


± 


0.20 


16.14 


± 


0.01 


-21.53 


-22.92 


2.25 


± 


0.08 


5.18 


1440+122 


0.162 


H 


16.93 


± 


0.12 


16.71 


± 


0.02 


-22.52 


-22.93 


3.90 


± 


0.25 


10.87 


1458+224 


0.235 


H 


15.78 


± 


0.08 


17.80 


± 


0.05 


-24.57 


-22.83 


3.20 


± 


0.80 


11.96 


1514-241 


0.049 


H 


14.48 


± 


0.12 


14.45 


± 


0.01 


-22.22 


-22.30 


3.70 


± 


0.10 


3.55 
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Table 4 — Continued 



Name z Type a nucleus 15 vciji host c ran nucl. d ran host e r e (arcscc)' r e (kpc) g 



1534+014 





.312 


H 


19. 


.08 


± 





.15 


18. 


16 


± 


0. 


02 


—21. 


.98 


-23. 


31 


2, 


00 


± 





10 


9. 


16 


1704+604 





280 


H 


21 


.08 


± 





.35 


18 


69 


± 


0. 


01 


-19. 


71 


-22 


45 


0, 


.85 


± 





.03 


3, 


61 


1728+502 





055 


H 


16. 


43 


± 





.10 


15 


49 


± 


0. 


02 


-20. 


53 


-21 


52 


3 


15 


± 





.05 


3, 


37 


1749+096 


0. 


320 


L 


16. 


.88 


± 





.05 


18 


.82 


± 


0. 


10 


-24. 


.24 


-22. 


73 


3. 


00 


± 





.80 


13 


.97 


1757+703 





.407 


H 


18. 


43 


± 





14 


19 


.58 


± 


0. 


25 


-23. 


31 


-22 


77 





.85 


± 





.50 


4, 


62 


1807+698 





.051 


L 


14. 


95 


± 





25 


13. 


.87 


± 


0. 


02 


-21. 


84 


-22. 


.97 


2, 


.10 


± 





.10 


2, 


09 


1853+671 





212 


H 


19. 


.48 


± 





.10 


18 


.19 


± 


0. 


01 


-20. 


.63 


-22 


17 


1, 


50 


± 





.08 


5. 


18 


1959+650 





.048 


H 


15. 


.38 


± 





.10 


14 


.92 


± 


0. 


02 


-21. 


.27 


-21. 


.78 


5, 


10 


± 





,10 


4. 


80 


2005-489 





071 


H 


12. 


.73 


± 





01 


14. 


52 


± 


0. 


01 


-24. 


81 


-23 


09 


5, 


65 


± 





,08 


7, 


65 


2007+777 





.342 


L 


18. 


.03 


± 





.10 


19 


.03 


± 


0. 


10 


-23. 


.26 


-22. 


73 


3, 


30 


± 





,90 


16 


.07 


2143+070 





237 


H 


18. 


.21 


± 





.11 


17 


.89 


± 


0. 


02 


-22. 


.16 


-22. 


76 


2. 


.10 


± 





,15 


7. 


90 


2200+420 





.069 


L 


13. 


.58 


± 





05 


15 


37 


± 


0. 


02 


-23. 


,89 


-22 


.17 


4, 


80 


± 





.40 


6, 


33 


2201+044 





.027 


L 


17. 


.18 


± 


0. 


05 


13. 


,74 


± 


0. 


01 


-18. 


.18 


-21. 


,64 


6, 


.78 


± 





,08 


3. 


68 


2254+074 





.190 


L 


16. 


.94 


± 





.12 


16 


.61 


± 


0. 


02 


-22. 


.90 


-23. 


,44 


4, 


90 


± 





.35 


15 


.54 


2326+174 





.213 


H 


17. 


63 


± 





.11 


17 


.56 


± 


0. 


03 


-22. 


.49 


-22 


,81 


1, 


.80 


± 





.15 


6, 


24 


2344+514 





.044 


H 


16. 


83 


± 





.05 


14. 


.01 


± 


0. 


01 


-19. 


62 


-22 


,49 


5, 


93 


± 





.02 


5, 


14 


2356-309 





165 


H 


17. 


.28 


± 





.13 


17 


21 


± 


0. 


02 


-22. 


.21 


-22. 


,47 


1, 


.85 


± 





.10 


5. 


23 



a SED type: H=HBL with F lk< , v /F 5G Hz > 5.5, and L=LBL, with F lkeV /F 5GHz < 5.5. 

b Apparcnt R-band magnitude of nucleus from best-fit PSF + de Vaucouleurs model. 

c Apparent R-band magnitude of host galaxy from best-fit PSF + dc Vaucouleurs model. 

d Absolute R-band magnitude of nucleus, K corrected assuming a = 1. 

° Absolute R-band magnitude of host galaxy, K corrected assuming early-type spectrum. 

f De Vaucouleurs effective radius in arcseconds. 

g De Vaucouleurs effective radius in kiloparsecs, assuming go = 0. 



Assumed 


Slope 


Intercept 


Zform 


(K-band mag/z) 


(K-band mag) 


2.8 


-1.37±0.30 


-24.97±0.11 


1.8 


-0.98±0.29 


-25.02±0.10 


1.3 


-0.65±0.28 


-25.04±0.10 



Table 5: Regression analysis results for absolute luminosity evolution based on the range of 
evolution model fits to the Hubble diagram. 
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Fig. 1. — Central regions of STIS F28x50LP images right, with contour plots left. 

Fig. 2. — Two-dimensional and three-dimensional images of the composite stellar PSF con- 
structed from archival STIS F28x50LP images of stars. 

Fig. 3. — Azimuthally averaged profiles of STIS images. The upper left and right panels show 
the best-fit de Vaucouleur + PSF models and exponential disk + PSF models respectively. 
Middle panels show zoomed-in profiles, inset with 1, 2 & 3a x 2 contours projected on m/i OSt -r e 
plane (for cases where the host galaxy was confidently detected.) Lower left panels show the 
best-fit dual-component models and lower right panels show the best-fit PSF-only models. 
For 1308+326, only the best PSF + exponential disk upper and lower left and PSF-only right 
profiles are shown, while for 2240-260 the best PSF-only profile is shown. 

Fig. 4. — BL Lac images with best-fit models subtracted — top: unsubtracted image; mid- 
dle: PSF subracted, normalized according to best-fit de Vaucouleurs model; bottom: best 
PSF + de Vaucouleurs model subtracted. For 1308+326 and 2240-260, the unsubtracted 
image (top) and best-fit PSF subtraction assuming no host galaxy (bottom) are shown. The 
arrow next to the object name indicates north and east. In the case of 0235+164, we chose 
a gray scale to bring out the detail in the surrounding structure, and so the marginal host 
galaxy is difficult to make out. 

Fig. 5. — Hubble diagram for BL Lac host galaxies, in F28x50LP magnitudes, corrected for 
galactic extinction. The large points are the high-z sample: STIS targets from this study. 
The small points are the low-z sample: WFPC2 snapshot targets with z < 0.6 from Urry et 
al. (2000). Also shown are the expected tracks of a non-evolving stellar population with an 
early-type spectrum (B — V = 0.96; dashed line), and a late-type spectrum (B — V = 0.57; 
dotted line). 

Fig. 6. — Hubble diagrams of high-z (large points) and low-z (small points) BL Lac host 
galaxies. Top shows the best-fit passively evolving model, with Zf orm = 1.8 (solid line), and 
the Zf orm = 5 passive model (dashed line). Middle shows the best-fit non-evolving model 
which also fits the low-z colour constraint of R — H > 2, with age = 5 Gyrs (solid line), 
as well as the age = 1 Gyrs (dashed line) and age = 0.02 Gyrs (dotted line) non-evolving 
models. This latter extremely young population is needed for a non-evolving population 
to fit the data. Bottom shows a dual population (solid line), in which 98% of the mass 
is evolving passively from Zf orm = 5, and 2% is undergoing active star formation with a 
constant age of 1 Gyr down to z — 0.6, after which it evolves passively. 
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Fig. 7. — Hubble diagrams of high-z (large points) and low-z (small points) BL Lac host 
galaxies with evolutionary models including dust and metallicity evolution. Top left shows 
a 0.8 Gyr non-evolving population with dust increasing linearly to an extreme optical depth 
of ry(0) = 8 between z = 0.8 to 0.4 (solid line), giving a reasonable fit to the data. The 
dashed line shows dust increasing to Ty(0) = 2, and the dotted line shows the model with 
no dust, both yielding poor fits. Top right shows a Zf orm = 5 passive population with dust 
increasing from r v (0) = to 10 by low redshift, producing a reasonable fit, while an increase 
to Ty(0) = 2 (dashed line) or a model with no dust (dotted line) produce poor fits. Bottom 
left shows the non-evolving population with metallicity increasing from 0.02 solar to 1 solar 
at z — 0.6. The solid line produces a reasonable fit with dust also increasing to Ty(0) = 5 by 
low redshifts. An increase to Ty(0) = 2 (dashed line) produces a marginal fit. The dotted line 
shows the model with no change in dust or metallicity. Bottom right is the same as bottom 
left, but with a Zf orm = 5 passively evolving population, which requires the same dust and 
metallicity evolution to fit the data. 



Fig. 8. — Absolute K-band magnitude versus redshift for high-z (large dots) and low-z (small 
dots) BL Lac host galaxies. These are K corrected assuming a redshift- dependent spectrum, 
determined using the evolution model which provided the best fit to the Hubble diagram 
(passive evolution with z/ orm = 1.8; see §4.3). Also plotted are the best-fit regression lines, 
derived using survival analysis to account for the upper limits. The solid line shows the best 
fit for Zf orm — 1-8 K corrections, while the other lines represent the 3a upper and lower limits 
to the slope, determined assuming Zf orm = 2.8 K corrections (dashed line) and Zf orm = 1.3 
K corrections (dotted line). 
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